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Understanding the nature of strong many-body correlations in condensed-matter 
finite systems is a fundamental task, which can be facilitated by devising analytical or 
semi-analytical many-body wave functions that approximate well the exact solutions 
and capture the essential physical properties of such systems. The description of strong 
correlations in the two-dimensional (2D) electron gas in a high magnetic field by the 
Laughlin wave function [Q represents a celebrated example of such a methodology. 

Here we construct semianalytical correlated wave functions for electrons in 2D 
single and molecular Quantum Dots (QD's), which are manmade devices created 
at semiconductor interfaces and are most often referred to as artificial atoms and 
molecules (since they contain a finite number of electrons). To this end, we employ 
a two-step approach consisting of symmetry breaking at the Hartree-Fock single- 
determinantal level and subsequent symmetry restoration via Projection Techniques 
(PT's) resulting in multi-determinantal wave functions. PT's have been introduced 
earlier in Quantum Chemistry |Q for the restoration of the total spin of a molecule and 
in Nuclear Physics |4| for the restoration of the total 3D angular momentum (space 
rotational symmetry) of deformed open-shell nuclei. Our application of these methods 
to circular single quantum dots (SQD's) requires the simultaneous restoration of both 
the spin and the angular-momentum symmetries. This requirement poses a more 
demanding challenge compared to the task of restoring a single broken symmetry, as 
is the case with the Quantum-Chemistry and Nuclear-Physics many-body problems. 

Since the approach of restoration of broken symmetries is largely untested 
and unknown in the framework of 2D QD's, we aim in this paper at focusing on 
methodological aspects. For simplicity and conceptual clarity, we construct below 
strongly correlated wave functions for the case of two interacting electrons in a SQD 
(artificial helium, He-QD) and in a lateral double quantum dot (artificial hydrogen 
quantum dot molecule, H2-QDM). Two strongly-correlated electrons (with or without 
an external magnetic field) may exhibit a complex behavior |^, |^ and can provide 
the foundation for understanding the properties of a larger number of interacting 
electrons. Furthermore, the recent growth of interest in quantum computing focussed 
attention on the potential ability to manipulate the ground-state entanglement of 
two-electron quantum-dot systems [Q. 

The two-body hamiltonian for two interacting electrons constrained to move on 
a plane is given by, 

Hiri) + H{r2) + eyKru , (1) 

where the last term is the Coulomb repulsion with k being the dielectric constant of 
the semiconductor material. H{r) is the single-particle hamiltonian for an electron in 
an external perpendicular magnetic field B and an appropriate potential confinement. 
For a QDM, the external confinement is given by a two-center-oscillator potential, and 
the single-particle hamiltonian is written as, 

H = T+ \m*i,l{x^ + y'i) + Kcck(y) + • s , (2) 

where y'^. ~ y ~ yk with fc = 1 for y < (left) and fc = 2 for y > (right), x 
denotes the coordinate perpendicular to the interdot axis (y). T — {p~ eA/c)^/2m*, 
with A = 0.5{—By,Bx,0), and the last term in Eq. (|^) is the Zeeman interaction 
with g* being the effective g factor, fig the Bohr magneton, and s the spin of an 
individual electron. The shapes described by H are two equal semi-circles connected 
by a smooth neck [Vneck(y)]- —yi = 2/2 > are the centers of these semi-circles, 
d = y2 — yi is the interdot distance, and m* is the effective electron mass. The case 
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of a SQD is obtained for d — 0, and in this case the confining potential is purely 
parabolic. For the smooth neck, we use Kcck(?/) = ^'m*ujQ[Cky'k + Dky'^]9{\y\ ~ IVkDi 
where 9{u) = for u > and 9{u) = 1 for u < 0. The constants Ck and D/. can be 
expressed via one parameter, as follows: Ck — {2 — 4e'')/?/fe and Dk — (1 — 3e'')/y^, 
where the barrier-control parameter = Vb/Vo is related to the actual (controlable) 
height of the bare barrier (Vb) between the two QD's, and Vq = to*WqJ/^/2. The 
single-particle levels of H are obtained by numerical diagonalization (for details see 
Ref. 0). 

In each case, in the first step of our procedure, the two-electron hamiltonian [Eq. 
(0)] is solved ||, |l^ in the (symmetry-breaking) spin-and-space unrestricted Hartree- 
Fock (sS-UHF) approximation; for comparison, the symmetry-adapted restricted 
Hartree-Fock (RHF) will also be considered. In all cases, we will use Hujo = 5 meV 
and TO* = 0.067TOe (this effective-mass value corresponds to GaAs). We will vary 
the dielectric constant k, however, in order to control the ratio of the strength of the 
Coulomb repulsion relative to the zero-point kinetic energy (see in particular the case 
of the SQD and figure 3 below). The variation of this ratio provides us with the ability 
to study the whole range of electron correlations, from the regime of weak correlations 
to that of strong correlations. 

We start with the case of a lateral H2-QDM in a magnetic field, where only 
the spin projection needs to be considered, since the confining potential lacks circular 
symmetry. The example we discuss here^ corresponds to the case of weak interelectron 
repulsion (the full choice of parameters is given in the caption of Fig. 1). The sS- 
UHF determinant (henceforth we drop the prefix sS in subscripts) which describes the 
"singlet" (see below) ground state of the H2-QDM is given by the expression, 

|<I>uhf(1,2)) = |?x(l)i)(2))/V2 . (3) 

In Eq. (||), we have used a compact notation for the Slater determinant with 
u(l) = u(ri)a(l) and v{2) = w(r2)/3(2), where u(r) and v{r) are the Is-like (left) 
and Is'-like (right) localized orbitals of the sS-UHF solution, and a and P denote 
the up and down spin functions, respectively. Such orbitals for the field-free case are 
displayed in Fig. 1(a) (left column). Similar localized orbitals (which are complex 
functions) can appear also in the B ^0 case [see Fig. 1(b), right column]. 

|*I'uhf(1, 2)) is an eigenstate of the z-projection of the total spin, S = S1+S2, with 
eigenvalue Sz = 0. However, it is not an eigenstate of the square, S^, of the total spin 
(thus the quotation marks in "singlet" above). From the determinant |<1'uhf(1, 2)), 
one can generate a singlet eigenstate of (with zero eigenvalue) by applying the 
projection operator Pq = 1 — von-, where the operator 13712 interchanges the spins of 
the two electrons. 

Thus the singlet state of the two localized electrons is given by the projected wave 
function, 

|*(1, 2)) = Po|<fuHF(l, 2)) (X \u{\)v(2)) - \u(\)v{i)) . (4) 

In contrast to the single-determinantal wave functions of the RHF and sS-UHF 
methods, the projected many-body wave function (^ is a linear superposition of two 
Slater determinants, and thus it is an entangled state representing a corrective (post- 
Hartree-Fock) step beyond the mean-field approximation. We notice that the spatial 
reflection symmetry is automatically restored along with the spin symmetry. 

X For a systematic study of the H2-QDM at different values of k [including n = 12.9 (GaAs)] and 
interdot barrier heights, demonstrating gate control of the entanglement of a pair of electrons, see 
Ref. 0. 
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Eq. (Q) has the form of a Heitler-London (HL) or valence bond ||Tl| wave function. 
However, unhke the original HL scheme which uses the orbitals <^L(r) and 4>r{v) 
of the separated (left and right) QD's, expression (^) employs the sS-UHF orbitals 
which are self-consistently optimized for any separation d, ineterdot barrier height 
Vb, and magnetic field B. Consequently, expression (|^) can be characterized as a 
Generalized Heitler-London (GHL) wave function. In the context of QD's, the simple 
HL approach has been proven very useful in demonstrating |Q that the H2-QDM 
under an external magnetic field can function as an elemental two-qubit logic gate. 
Our more accurate GHL approach has the potential of greatly improving the quality 
of similar investigations regarding the implementation of quantum computing using 
solid-state nanodevices. 

We further notice that our GHL method belongs to a class of projection techniques 
known as variation before projection (VBP), unlike the familiar in Quantum Chemistry 
Generalized Valence Bond method (l2|, or the Spin-Coupled Valence Bond method 
[I3[, which employ a variation after projection (VAP) (see e.g. ch. 11.4.2 of Ref. 
[p). In the context of QD's, a generalization of our GHL approach along the VAP 
technique may provide even more accurate results. This, however, is left for future 
work, including the development of the pertinent computer codes. 

The energy of a projected state is given Q in general by the formula, 

£^PROJ = ($uhf|HO|$uhf>/($uhf|0|$uhf) , (5) 

where Ti. is the many-body hamiltonian and O is a general projection operator with 
the property = O. 

Using the spin-projection operator Pq in place of O, we find for the total energy, 
-^GHL' '^^ singlet GHL state, 

+ Juv + Kuv] , (6) 

where H^m etc., are the matrix elements of the single-particle hamiltonian H in Eq. 
(^, and J and K are the direct and exchange matrix elements of e'' j Kr\2- Suv is the 
overlap integral of the u{r) and v{r) orbitals, 

Suv= j u*{Y)v{v)dv , (7) 

and the normalization parameter is given by 

AC' = 1/(1 + Su.S^u) ■ (8) 

At zero magnetic field, the electron orbitals are real functions and Eqs. (||) — 
(^ reduce to a form familiar from the Generalized Valence Bond theory of Quantum 
Chemistry |l|l. 

For the triplet state with = ±1, the projected wave function coincides with the 
original HE determinant, so that the corresponding energies in all three approximation 
levels are equal, i.e., -Eqhl = ^rhf = ^uhf- 

In Fig. 1(c), we display the singlet-triplet energy gap, Ae = E^ — E^, of the 
H2-QDM as a function of the magnetic field B. For all three levels of approximation, 
Ae starts from a negative minimum (singlet ground state) and after crossing the zero 
value it remains positive (triplet ground state). However, after crossing the zero line, 
the RHE curve incorrectly continues to rise sharply and moves outside of the plotted 
range. After reaching a broad maximum, the positive Ae branches of both the sS-UHF 
and GHL curves approach zero for sufficiently large B, a behavior which indicates that 
the H2-QDM dissociates with B. This molecular dissociation can be further seen in 
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FIGURE 1 
SEPARATE GIF 



Figure 1. H2-QDM: The two occupied orbitals (modulus square) of the 
symmetry broken "singlet" sS-UHF solution for (a) B = and (b) B = 9 T. 
(c) The singlet-triplet energy difference as a function of an external perpendicular 
magnetic field and for three successive levels of approximation, i.e., the RHF (top 
solid curve), the sS-UHF (dashed curve), and the GHL (bottom solid curve). 
The choice of parameters is: parabolic confinement of each dot Hluo = 5 meV, 
interdot distance d = 30 nm, interdot barrier height Vj = 3.71 meV, effective 
mass m* = 0.067me, and dielectric constant k = 45. The effective Zeeman 
coefficient was chosen g* = 0; for a small nonvanishing value, like g* = —0.44, the 
Zeemann splitting does not alter the orbital densities and can be simply added to 
Ae. Up and down arrows denote spins. Distances are in nm and the densities in 
10-* nm-2. 



the orbitals themselves: at i? = 9 T [Fig. 1(b)] the orbitals are well localized on 
the individual dots, while at B — [Fig. 1(a)] they extend over the entire QDM. 
Note that, in addition to having the proper symmetry, GHL is energetically the best 
approximation. 

We further note that systematic explorations for determining the border of HF 
instability (i.e., the appearance of broken symmetry solutions with lower energy) in 
the case of the H2-QDM are not available. At i? = 0, this border depends on all four 
parameters hujo, k, d, and Vb- For Hloq = 5 meV, d = 30 nm, and k = 45, it was found 
[ pO[ that this border can be crossed by reducing the interdot barrier Vb to zero. For 
the same values of B, Tiloq and d, but for a different k — 12.9 (stronger interelectron 
repulsion, as is the case for GaAs), however, all values of the interdot barrier {Vb > 0) 
represent cases that lie well within the region of instability. 

Next we address the case of two electrons at B=0 in a single QD with a parabolic 
confinement. As we have shown in earlier publications jH |9|, the sS-UHF solution 
does not preserve the rotational symmetry for Rw > 1.0 (for Rw < 1.0, the sS- 
UHF solution collapses onto the RHF one). The Wigner parameter Rw expresses 
the ratio between the interelectron repulsion Q and the zero-point kinetic energy 
K; it is customary to take Q — e^/ulo and K = hujQ, where Iq = [h/m*LOQ)^/'^ is 
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the spatial extent of an electron in the lowest state of the parabolic confinement. 
This gives Rw oc l/{Ky/LJ^), showing that it can be varied through the choice of 
materials (i.e., k) and/or the strength of the confinement {ujq). For Rw > 1-0 and 
for Sz = 0, the sS-UHF yields |^ an infinite manifold of azimuthally degenerate 
ground-state configurations consisting of two localized antipodal orbitals, suggesting 
formation of an electron molecule ||, ^ (also referred to as Wigner molecule). For 
K = 8 {Rw = 2.39) and for a specific azimuthal direction (i.e., 7 = 0), the two 
localized orbitals, u{r) and v{r) with up and down spin, respectively, are depicted in 
Fig. 2(a). As was the case with the double-dot example studied above, the Sz — 
sS-UHF determinant |$uhf(1,2)) for 7 = is given by Eq. (||) and thus, in addition 
to the rotational symmetry, it does not preserve the total-spin symmetry. For the 
singlet state, one can restore both symmetries successively, namely, one can generate 
appropriate projected wave functions, 

|vl^/(l,2)) =0|$uhf(1,2)) , (9a) 

by applying the product operator ,| 

O = ViPo , (96) 

where the spin-projection operator Pq produces a two-determinant singlet GHL wave 
function as previously explained. The angular-momentum projection operator Vi 
yields multideterminantal wave functions having good total angular momentum /. 
This latter operator produces a linear superposition of an infinite number of azimuthal 
GHL wave functions and it is given|| by 

2ttVi= / d7exp[-i7(L-/)] , (10) 
Jo 

where L — li ^- 12 is the operator for the total angular momentum. In the following 
we focus on the ground state^ of the system with / = 0. 

It is instructive to examine the transformations of the ground-state electron 
densities (ED's) resulting from the successive approximations, RHF, sS-UHF, spin 
projection (SP), and combined spin and angular momentum projection (S&AMP). 
For K = 8, such ED's are shown in Fig. 2. Since the exact solution for two electrons 
in a parabolic confinement is available we also plot the corresponding ED for 
the exact ground state in Fig. 2(f). The ED's of the initial RHF [Fig. 2(b)] and 
the final S&AMP [Fig. 2(e)] approximations are circularly symmetric, while those of 
the two intermediate approximations, i.e., the sS-UHF and SP, do break the circular 
symmetry. This behavior graphically illustrates the meaning of the term restoration 
of symmetry and the interpretation that the sS-UHF broken-symmetry solution refers 
to the intrinsic (rotating) frame of reference of the electron molecule. We notice that 
the S&AMP electron density exhibits a characteristic flattening at the top in contrast 
to the more Gaussian-type RHF one. Thus, although not identical (the exact ED is 
slightly flatter at the top), the S&AMP ED closely resembles the exact one displayed 
in Fig. 2(f). Further, we remark that the SP electron density exhibits in the middle 

§ In the case of a triplet state, the sS-UHF determinant with Sz = ±1 preserves the total spin, and 
thus the application of Vj alone is Mifficient. 

II The corresponding expression pi for the 3D angular-momentum projection uses the Wigner 
functions ©^^^(Q) (see also ch. 11.4.6 in Ref. [Q). 

1 The family of projected wave functions (9) describes all the lowest-energy (yrast band fel) states 
with good angular momentum / = 2, 4, in addition to the ground state (7 = 0). The yrast-band 
states with odd values, 7=1,3, 5, are generated via a projection of the triplet state. 
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FIGURE 2 
SEPARATE GIF 



Figure 2. Hc-QD at B = 0: (a) The two occupied orbitals (modulus aquare) of 
the symmetry-broken "singlet" sS-UHF solution. The rest of the frames display 
the electron densities at successive approximation levels, i.e., (b) RHF, (c) sS- 
UHF, (d) SP alone, (e) S&AMP, and (f) exact. See the text for a description of 
the acronyms. The choice of parameters is: dielectric constant k = 8, parabolic 
confinemenet huio = 5 meV, and effective mass m* = 0.067me. Distances are in 
nm and the densities in 10~* nm~^. 



a shallower depression than the corresponding sS-UHF one, in keeping with the fact 
that the covalent bonding increases the probability of finding the electrons between 
the individual dots. 

Using the projection operator (^) in the general Eq. (|^), we obtain for the energy 
EsSzAMP of the fully projected ground state of the two-electron single QD, 

£^S&AMP = j h{l)dl I J n{'j)d'y , (11) 

with 

h{'j) = HusSyt + HutSys + HytSuS + H V S S ut + VuV St + Vuvts , (12) 

and 

n{"f) = SusSvt + SutSys , (13) 

s(r) and t{r) being the u(r) and v{r) sS-UHF orbitals rotated by an angle 7, 
respectively. Vuvst and Vuvts are two-body matrix elements of the Coulomb 
repulsion.]^ Observe that i?s&AMP represents a rotational average over all the possible 
orientations of the electron molecule and that /i(0)/n(0) coincides with the expression 
(^ for the energy i?QHL of the covalently bonded H2-QDM. 

+ Vuvst = (e^/k) / dri J dr2U*(ri)u*(r2)(l/ri2)s(ri)t(r2). In Eq. J„„ = Vuvuv and Kuv = 
Vuvvu. 
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Figure 3. He-QD at S = 0: Difference between the ground-state energies 
of various approximations and the exact one, plotted vs. the dielectric constant 
ft. From top to bottom, the curves correspond to the RHF, the sS-UHF, the 
SP alone, and the combined S&AMP approximation. The curve labeled RHF 
gives the correlation energy. The choice of parameters is: parabolic confinement 
hujo = 5 meV, and effective mass m* = 0.067me. 



To further examine how closely expression (9) describes the ground state of the 
He-QD, we display in Fig. 3 the energy deviations of the four approximation steps (i.e., 
RHF, sS-UHF, SP, and S&AMP) from the exact (EX) [§ ground state, as a function 
of the dielectric constant k. The variation of k produces a variation in the strength of 
the interelectron repulsion, with a larger repulsion (i.e., smaller k) corresponding to 
stronger electron correlations. 

The correlation energy, i?cor = ^-rhf — ^'EX, is defined as the difference between 
the RHF ground-state energy and the exact one, and thus it coincides with the top 
curve in Fig. 3 denoted as RHF. The ordering (from top to bottom) of the curves in 
Fig. 3 reflects the fact that each subsequent approximation step captures successively 
a larger portion of the total correlations, and thus lowers the total energy. From Fig. 
3, one sees that there are two correlation regimes: In the regime of weak correlations, 
the three lower curves collapse onto the RHF one; naturally, this regime corresponds 
to the normal Fermi liquid and extends to k — s- oo. The onset of strong correlations 
and of symmetry-broken HF solutions occurs at k = 16.5 {Rw — 1-16) and extends 
to /t = 0. 

Note that the percentage, Ecoy/Eex, of correlations with respect to the exact 
energy is 6.8% and 49.5% for k = 16.50 and k — 1, respectively. These values 
are much higher than the values encountered in natural atoms; for this reason, 2D 
QD's are referred to as strongly correlated systems. Thus it is not surprising that 
the RHF error in Fig. 3 grows exponentially for stronger correlations (smaller k or 
larger Rw)- However, it is remarkable that the S&AMP approximation (where all 
symmetries have been restored) converges to the exact result for smaller k. Note that 
the exact ground-state energy is 19.80 meV for k 8 {Rw = 2.39) and 51.83 meV 
for K = 1 {Rw = 19.09), resulting for the S&AMP in a relative error of 4% and 
0.7%, respectively-!^ These results are in agreement with the experience from Nuclear 
Physics, where it has been found [Q that the VBP yields reliable results in the case 

* For K = 8, the fraction of correlation energy captured by the successive approximations is: sS-UHF 
44.4%, SP 57.2%, and S&AMP 73.1%; the corresponding values for k = 1 are: 94.5%, 94.5%, and 
98.5%. 
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of strong symmetry breaking (strongly deformed nuclei). 

Interestingly, the SP curve collapses onto the sS-UHF one for k < 3.2. This 
behavior defines Q an intermediate regime|| between the Fermi-liquid and the regime 
of strongly crystallized Wigner molecules (k < 3.2). In the latter regime, the overlap 
(Suv) between the antipodally localized electron orbitals is very small, such that the 
sS-UHF energy is not effectively lowered by the spin projection (SP). 

We have demonstrated that the sS-UHF method in conjunction with the 
companion step of restoration of the spin and space symmetries via Projection 
Techniques (when such symmetries are broken) can provide appropriate semianalytical 
wave functions for an accurate description of strongly correlated electrons in artificial 
atoms and molecules. 

This research is supported by the US D.O.E. (Grant No. FG05-86ER-45234). 
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